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Abstract. We consider two variants of a quantum-statistical generalization 
of the Cramer-Rao inequality that establishes an invariant lower bound on the 
mean square error of a generalized quantum measurement. The proposed com- 
plex variant of this inequality leads to a precise formulation of a generalized 
uncertainty principle for arbitrary states, in contrast to Helstrom's variant 
in which these relations are obtained only for pure states. A notion of canon- 
ical states is introduced and the lower mean square error bound is found for 
estimating of the parameters of canonical states, in particular, the canonical 
parameters of a Lie group. It is shown that these bounds are globally attain- 
able only for canonical states for which there exist efficient measurements or 
quasimeasurements. 



The development in recent years of the theory of generalized quantum mea- 
surements (see the review and the literature cited there) has made it possible 
to introduce the concept of a quasimeasurement of incompatible observables de- 
scribed by noncommuting operators and, using this, to solve a number of problems 
of the quantum theory of information and communication [3J 0] |5J |BJ |S| , give for 
pure states a precise formulation of a generalised Heisenberg uncertainty principle 
for quantities such as, for example, the time and energy, or phase and number of 
quanta [Hj, and to define precisely what is a measurement of the time and phase 
in quantum mechanics 0- In accordance with this theory, every quantum 
measurement in this generalised sense is described by a positive resolution of the 
identity operator 1 on the Hilbert space H of state- vectors \ip) of the observed 
quantum system: 



Here II (•) is an additive mapping (measure) on the Borel algebra 05 (X) of a mea- 
surable space X 3 x into the set of Hermitian-positive (i.e. nonnegative-dcfmitc 
Hermitian) operators in TL. Such normalized positive measure II will be called 
quantum probability measure ( QPM), or simply quasimeasurement. If g is a quan- 
tum state density operator, the probability Pr (B) of an event x E B in such a 
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measurement is evaluated in accordance with the formula 



Pr (B) = TrpIT (B) , Be<8 (X) 



where Tr denotes the usual trace in TL. If the quantum measure II in is 
orthogonal, II (A) II (C) = for every A n C = 0, then it is a projerctor-valued 
measure. The generalised measurement in this case with X = M. n reduces to an 
ordinary measurement of the commuting self-adjoint operators 



For the nonorthogonal QPM, there is no one-to-one correspondence between (|l.lf) 
and (|1.2fl . The corresponding generalized measurements, which are called hence- 
forth approximate measurements of the operators l|1.2fl . are not described as the 
measurements of these Hermitian operators even if they commute, though fre- 
quently they can be described uniquely by a single non-Hcrmitian (non-normal) 
operator (see Section or, more generally, by a family of noncommuting Hermit- 
ian operators. 

There is an intimate connection between the concept of a quasimeasurement as 
approximate measurement and the concept of an indirect quantum measurement 
(an indirect measurement is an ordinary measurement in an extended quantum 
system that includes the original system as a part 3 ). This connection is a sim- 
ple consequence of the Naimark's well-known theorem on the existence for every 
nonorthogonal QPM of an orthogonal one in an extended Hilbert space that com- 
presses to the original QPM on the subspace TL. 

One of the results of this paper is to show how the concept of a generalized 
measurement enables us to formulate precisely a generalised Heisenberg uncertainty 
principle for quantities such as the time and energy, phase and number of quanta, 
angle of rotation and angular momentum as a consequence of a quantum Cramer- 
Rao type inequality for the arbitrary states. The first members in each of these pairs 
- the time, phase, and angle - cannot, as is well known, be described by Hermitian 
operators in TL, though their measurement can be described as a statistical estimate 
of the corresponding parameters of quantum states. As Helstrom has shown in |S] 
by means of the symmetric quantum Cramer-Rao inequality which he introduced 
in the variances of the results of any measurements to obtain such an estimate 
for pure states cannot be lower than a certain level that is inversely proportional to 
the variances of the generators of the unitary representations of the corresponding 
translation groups (i.e., the operators of the energy, number of quanta, or angular 
momentum). For example, if a pure state of a harmonic oscillator is known up to 
the oscillator phase, its state- vector is unitarily equivalent to a fixed vector I^q) G TL 
and can be described by the family 



where n, the operator of the number of quanta, is generator of the representation 
e i8n/h Q f g r0U p f phase translations. If a QPM n determines the probabilities 



(1.2) 




= e' 
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it defines an unbiased estimate of the unknown value of the phase 9 as the mea- 
surement result A. The corresponding quasimeasurement, described by the induced 
QPM IT (dA) on (— 71", tt], is an approximate measurement of the "phase operator" 
q = J ALT (dA) .The mean quadratic error of the measurement approximation for 
such q is given by the quantum expectation {ipg\a 2 \^g) of the positive operator 

a 2 = J (X-q)U^ (dA) (X-q), 

and the total variance describing the estimation accuracy of 9, 

(A — Q) 2 , 

is the sum of this and mean square distance of the operator q and 91: 

Re = ^e\(t 2 + (<i-0) 2 )\^e)- 

The quantum Cramer-Rao inequality proves in this case that the second variance 
cannot be below the level h 2 /4 (^(n — ng) 2 ^ , and thus Rg > h 2 /AGg, where 

Ge = (ipe\ (" ~ n e) 2 \i>$) = ((« - n °) 2 ) ' 

is Fisher information as the variance of n with ng = (^>g\n \ipg) = no- This is 
Helstrom's precise formulation of the generalized Heisenberg's uncertainty principle 
for the conjugate quantities 9 and h, the first of which is described by a generalised 
measurement satisfying the unbiased condition 1)1. 3)1 . 

In Section [3 we give the invariant formulation (|2.4(l of the Helstrom's Cramer- 
Rao inequality, and we also consider another generalization (|2.7|) of this inequality, 
which in contrast to Helstrom's can be naturally adapted to a complex situation 
and enables one to obtain straightforward a multidimensional generalization of the 
uncertainty relations 1)3.5)1 for not only pure but also mixed states. We also obtain 
the noncommutative generalization 1)3. 12J) of these relations for the generators and 
canonical parameters of unitary representations of an arbitrary Lie group. These 
generalizations are intimately related to the canonical families of states described 
in Section |3J whose particular role is disclosed in Section in which it is shown 
that if the lower bounds for the mean square errors of a measurement are to be 
attainable, it is necessary and sufficient that the corresponding density operators 
have the canonical form 1)3.1(1 . 

2. Invariant Bounds of the Cramer-Rao Type in quantum Statistics 

1. Let {q$,$ £ M} be a family of density operators g$ in TL that describe 
the statistical state of a quantum system as a smooth function of unknown real 
parameters t9 = 0& , . . . ,$ m ) in a given manifold M C R m . Every simultaneous 
measurement of these parameters can be described in TL by QPM n which defines 
a row-vector random variable A 6 R m with probability distribution Pr(dA|"i?) = 
TrgjjII (dA) known up to d. The mean quadratic errors of the measurement are 
determined by the components 

(X 1 - (A fe - 



(1.3) 



Re := / (A - 9) Pr [dx\9) = 



Rf = M,, 



4 



V P BELAVKIN 



of the covariance matrix R$ = [R l Jf] by means of expressions trC T i? = CikR %k given 
by a Hcrmitian-positive matrix C = [cik] of the quadratic cost form 

c(A,0) = (V-^Cife (A fe -tf fc ) 

which plays the role of a metric tensor. Here and in what follows the Einstein 
summation convention is assumed: 

CikKtf = ^ ^ CikR^ '■ 

i k 

In what follows we shall consider only those measurements that satisfy the unbiased 
conditions M,j (A 1 ) = d l . under which the matrix is the covariance matrix of 
the estimates 0j, and the mean square error for fixed R$ takes a minimal value. 

Helstrom established £Q for the covariance matrix R$ of such measurements 
a lower bound by using the concept of operators gi of symmetrized logarithmic 
derivatives of the function q# with respect to i9\ He defined these gi by means of 
the equations 

Q 

(2-1) giQ# + QftQi = 2diQ € , d t := —. 

art 

As in the classical case this bound is determined by the matrix G,j = [Gik (i?)] 
of the covariances of the logarithmic derivatives §i = fji (1?) of Eqs. (|2.ip . defined 
in the symmetrized form as 

(2-2) G ik (0) = i (& (0) 3 fc (0) + & (1?) & 

(Note that due to Trd*^ = 

(9i (<% := Tr^.g 4 (0) = 
for all z?). The corresponding inequality has the form 
(2.3) R& > G^ 1 , e M, 

and is understood in the sense of nonnegative definiteness of the matrix [Rff — Gjf] , 
where G l $ are the components of the inverse matrix G^ 1 : G^Gjk = S k . The in- 
equality (|2.3() establishing an uncertainty relation between the variances of estima- 
tion and the variances of the corresponding logarithmic derivatives, is a quantum 
analog of the Cramer- Rao inequality . The matrix G# which we call symmetric 
quantum Fisher information, or more fair, Helstrom information, is one of possible 
generalizations of classical Fisher information. It plays the role of a metric tensor 
that locally defines the geodesic distance 

s G (0, + dtf) = (G ik (0) d?f di? fc ) 1/2 

in the parameter space M C K m ; this is analogous to Fisher information distance 
in classical statistics. 

2. Further, we shall consider a slightly more general situation in which the 
state parameters are not the measured parameters $ l but local coordinates a = 
(a 1 , . . . , a") of a smooth manifold 6 parametrizing the unknown g = g (a). The 
measured parameters are assumed to be known smooth functions 0j (a) of the 
unknown parameters a k . The corresponding generalization of Helstrom's inequality 
(12. 3|) is a lower bound for the matrix R (a) = [Rik (a)] of the covariances 

R lk (a) = M [(A, - 00 (A fc - \a] 
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of the estimates Xi in the form 

(2.4) R(a) > D (a) G (a) -1 D (a) T . 

that is invariant under the choice of the state coordinates a — (a fc ) . Here D (a) 
is the matrix [D^ (a)) of the partial derivatives Di k (a) = d'di/da k , D T = [d^k], 
and G (a) = [Gi k (a)] is the symmetric quantum Fisher information corresponding 
to the coordinates a, that is a matrix of the covariances 

Gm (a) = Tr [g k (a) ■ g t (a) g (a)} , g k ■ g t = - (g k gi + gifjk) , 



2 



of the Helstrom's logarithmic derivatives 

d 

Q (") 9k + 9kQ (") = Q (a) 

with respect to the coordinates or. 

The inequality (|2.4|l reduces to the classical Cramer-Rao inequality only when 
the family {g(a)} is commutative. For noncommutative families, one can have 
other quantum generalizations 0] , [S] of the Cramer-Rao inequality based on other 
definitions of the logarithmic derivatives; these lead to other lower bounds for R 
that may differ from Helstrom's invariant bound DG _1 D T . Moreover, in the non- 
commutative case it makes sense to consider also the complex-valued parameters 
as any quantum state has the natural complex coordinatization g — aa* /Tra*a 
in terms of the complex Hilbert-Schmidt operators a with is the adjoint operators 
a*as their complex conjugated. In the case of complex parameters £ C repre- 
sented by analytic functions di (a, a*) the particular importance is acquired by the 
following invariant generalization of the Cramer-Rao inequality based on the right 
and left logarithmic derivatives which were proposed independently by the author 
[Hj and Yuen and Lax 0. 

3. Suppose the parameters a k are given in pairs (^J k ,0 k ^j G R 2 which are 

complexified as }f^ k + i6 k = f3 k . Such parameters a £ K 2 ™, considered as complex 

n-columns, will often be denoted as — {p k ^j € C™, with 7 = (3 + (3 S M™ and 

6 = Im/3 e R n . The partial derivatives d k = d/df3 k , d k = d/dp k are defined by 
means of the partial derivatives d/dj l , d/dO % in the usual manner: 

_d_ _(J>_ _i_<M d f d ; d 
dp k ~ \dj k +1 d6 k ) ' d~p k ~ \ d l k ~ 1 d9 k 

such that d k (3 l = S l k = d k (3 l and d k fi l = = B k j3 l . 

The estimated parameters Di, ! = l,...,mas functions of complex a, a can still 
be real functions of 7 and 9. They are not assumed to be analytic with respect 
to a, but differentiable independently with respect to a and a (e.g. given by bi- 
analytic functions di (a, a') at a' — a). We define the non-Hermitian right and left 
logarithmic derivatives of the density operator g(a,a) by the relations 

(2.5) gh k = -^, ^ g= ~H^> k=l,...,n. 

The operators h k = h k (a, a) of the right derivatives with respect to a k are Her- 
mitian conjugate at each a to the operators h% = h k (a, a)* of the left derivatives 
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with respect to a k , and they both have zero expectations 

Tr/ifc (a, a) g(a,a) = = Trft,£. (a, a) g (a, a) . 

The corresponding quantum Fisher information is given by the matrix H = [Hki] 
of covariances 

(2.6) H ki (a, a) = Tr h k (a, a) hi (a, a)* g (a, a) . 

Obviously this matrix is Hermitian-positive, and under the assumption of its non- 
degeneracy it defines a positive-definite metric 

dsn = H k ida k da l 

in some complex domain O C C™ of the unknowns a £ O. 

4. Suppose a simultaneous measurement of the parameters i?j is described by 
a QPM Ilonl that determines the estimates Xi of $j as complex-valued random 
variables of x G X with respect to the distribution 

Pr [dA | a, a] = Trll (dA) p (a, a) 

parametrized by a. 

The mean quadratic errors of the measurement are determined by the matrix 
R (a, a) = [Rij (a, a)] of covariances 

Rij (a, a) = M [(A, - ^) (A, - 0,-) |a, a] 

which can be written as the sum i£y = (o^-) + Qij of two kind errors. The first 
one is given by the Hermitian-positive matrix of the elements 

(a}j) (a, a) = Tr [a^g(a, a)] 

as the quantum expectation of the covariance operators 



-2 
<?ij 



(A i -g i )H(dA) (Xi-qj) 



for the quantum estimates 



& = / A,n (dA) 

The second forms the mean quadratic error matrix Q = [Qij] 

Qij (a, a) = ((& - tf L ) (q* - dj)) (a, a) 

for the the operators qi "estimating" the parameters t?£. 

Assuming the convergence of the integral defining qi, the unbiasness condition 

M [di\a,a] := J Aj Pr (dA | a, a) = i? 4 (a, a) , 

for the estimates A^ can be written in the form of quantum unbiasness 

(qi) (a, a) = Tr^g (a, a) = i? t (a, a) . 

Under this assumption the matrix Q (a, a) — [Qij (a, a)] is the covariance matrix 
of the operators qi, and as it is shown in the Appendix, it has lower bound Q > 
DH _1 Dt, and therefore 

(2.7) R (a, a) > D (a, a) H (a, a) -1 D (a, a) f , 

where D = D (a, a), as in 1)2. 4|l . is the matrix [dik] of the derivatives d$i/da k , and 
= [9fc??iJ is the Hermitian adjoint matrix. 
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As we shall see, even in the real case # = if, the bound may lead to a 

lower bound that differs from Helstrom's bound 1)2.4(1 . We shall say that 1)2.7)1 is the 
right lower bound. Besides this bound, we can consider other bounds, for example, 
the "left" bound, which is based on the left logarithmic derivatives with respect to 
a. All these bounds are proved in the same way as 1)2.7)) see the Appendix. Note 
that the right bound in 1)2. 7[) is invariant under the change of variables (a fc ) i— > 
by replacing the derivatives with respect to a k by derivatives with respect to the 
new variables di = i?j (at, a) only under the analyticity condition ddi/da k — of 
the transforming functions t?j (a, a) = i?j (a) and the condition of nondegeneracy 
of the matrix of the derivatives ddi/da k . Therefore, the inequality 1)2.7)1 and its 
noninvariant form R > H _1 arc not equivalent unless not only the nondegeneracy 
of the matrix D but also the analyticity condition d'di/da k — (i.e., the condition 
that the functions i?j (a, a) are independent of a) hold. 

3. Canonical States and Uncertainty Relations 

In classical mathematical statistics, a particular role is played by canonical, or 
exponential, families of probability distributions, for which the Cramer-Rao bound 
is attainable for a special choice of the parameters In Sectional we shall show 
that in quantum statistics an analogous role is played by the density operators of 
the form 

(3.1) g(P,P) =x(My 1 ^ hAl Qo^ kAh , 

where Xk, k = 1, . . . , n are linearly independent operators in 7i, which may be non- 
Hcrmitian: x\ ^ Xk, and even need not commute with the conjugates x%x\ ^ x" k Xi. 
We shall assume that the generating function 

(3.2) X {M)=Trg e^e^, 

of the moments of these operators in the state g = g is defined in an open neigh- 
borhood of the origin /3 = of the complex space C™ with finite first and second 
moments 

^XU=0 = = J-yI^o, ( Xi xt) = ^-Ajn^ 

(the operators Xk in ()3.1)l can always be chosen to have zero expectations (xi} — 
9jlnx|^ = o = in the state g ). We shall call that the family of density opera- 
tors ()3.1)l canonical, with the parameters (3 k canonically conjugate to the quantum 
variables ik ■ In contrast to the classical case, even for selfadjoint Xk one can mean- 
ingfully consider complex values of the conjugate parameters j3 k . 

Particular interest attaches to the case, which does not have a classical analog, of 
the canonical states 1)3.1)1 when f} k are imaginary, (3 k — i8 k , and Xk are selfadjoint, 
Xk = Sfc = x\. The parameters 6k — h9 k (h is Planck's constant) then take the 
dimension and meaning of the classical variables which are dynamically conjugate 
to their shift generators §k- For example, :9 : , if s is the Hamiltonian, then 9 is the 
time, if s is the momentum then 9 is the position, and if s is the number of quanta, or 
angular momentum, then 9 is the phase, or polar coordinate. For (3 k = iftr 1 9k = f3g 
the canonical states g g — p (Pg,f3 g ) ()3.1)) become unitary equivalent 

(3.3) g g = e ieks */ h g Q e- iek '»/ h 

to the state g = g (0, 0) corresponding to the zero value 9 = 0. 
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Now we shall see that the inequality l|2.7|) with 9k = Mm/3 applied to the 
canonical family (|3.3|) with commuting §k — s£ immediately provides the precise 
formulation of a generalised Heisenberg uncertainty principle for an unbiased es- 
timation of 0. In this the right and left logarithmic derivatives with respect to 
/3 and /3 for the family i|3.1[l are equal to the symmetric logarithmic derivatives 
9k (7) = h- Mfc (7) with respect to 7 = f3 + 0. 

d d 
h k = s k - -^-X (0+0)= 9k = s k - tt^-X (0 + 0) = K- 

This implies that the Fisher informations H ((J, 0) and G (x) coincide with the 
covariance matrix S (/3 + 0) of the commutative family Sk given at the state ijQJl 

by 

(3.4) S lk =x((3 + py 1 (fa - fid eiP+PY*' (h - Mfe )) Q . 

While the complex quantum Cramer Rao bound (|2.7|l for the unbiased estimation of 

real parameters i?i (7) with 7 = /3 + (3 coincides in this case with the Helstrom's 
invariant bound (|2.4|) . it also gives immediately the uncertainty relation 

(3.5) R e >~fi 2 So\ S = S(0) 

for the unbiased estimation of 9i = H9 l based on the imaginary parts Q % = Im for 
the canonical coordinates (3 with the fixed 7 = 0. 

Indeed, setting Oi = Mm/3 1 such that d0i/d(3 k = Mjfe/2i, we obtain from (|2.7() 
the generalised Heisenberg uncertainty relation in the form 

R > \h 2 s~ 1 ( P - + 0). 

Here R is the mean quadratic error matrix 

R = M [(A, - Mm/3 1 ) {\ 3 - Mm/3 J ) \(3,0\ 

of unbiased estimates r ) i and S is the matrix of the covariances .4|> defining the 
uncertainty relation Ij3.5[l at7 = /3 + /3 = 0. 

The uncertainty relation (|3.5|l acquires the following matrix meaning: The covari- 
ance matrix Rg — R (/3 e , f3 e ) of the unbiased estimates for the canonical parameters 
6i of the translation group represented in Ti by the unitary transformations (|3.3|) 
with the selfadjoint generators §k is in the canonical uncertainty relation with the 
covariance matrix of these generators, 

Sik (0) := TrggSiSk = TrQ SiS k = S ik (0) , 

in the initial (and any other transformed) state g = g (0). 

This uncertainly relation holds for all commuting Hcrmitian operators Sj, not 
only for those like momenta which have dynamically conjugate observables §i. Hel- 
strom derived this generalized uncertainty relation l|3.5|l in one dimensional version 
from his bound for the particular case of pure states g = \ip ) (ip \ |SJ. However 
for this purpose the symmetric inequality l|2.3() is inappropriate, and this is why his 
derivation involved so complicated matrix elements calculations. 

The uncertainty relations naturally correspond to not symmetric but antisym- 
metric logarithmic derivatives, defined as the Hermitian solutions f>k = Pk (0) of 
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the von Neumann equations 

r ~ i ~ K d 

[Qe,Pk : = BePk - PkQe = -^~Be 

For the canonical family (|3.3[) we have the solutions pk (9) = Sfc which are uniquely 
defined by the condition Trg g pk (9) = 0. Assuming that the solutions pk (9) exist 
for an arbitrary parametric family g g , one can derive the generalized uncertainty 
relation for the covariance matrix R g of the unbiased estimates 



Me[0i] = (qi)e = e h & = J % n(d9) , 



in terms of the new quantum Fisher information matrix S g = [Sik (0)] given by the 
symmetric covariances 

S ik (0) = Tr& (9)- Pk (9) g g . 
It simply follows form of matrix inequality 

Re > Qe > ^S e \ 
where Qe = [Qik (9)] is the matrix of covariances 

Qik (9) = ((& - Oi) (q k - 9 k )) e 

with R e - Qe = j [(A, - &) U (dA) (A fe - q k )] > 0. 

Indeed, due to the unbiasness (q) g — 9 we have mean canonical commutation 
relations 

d 

(m,Pk (0)]) e = Trqi [Pk {9) , g e ] = \Qi)e = ihS ^- 

From this and ([q,p]) g — 2lm(qp) g , where q = q — 9, p = p — [A, we derive Q > 
fi 2 /S' _1 /4 by Schwarz inequality and \(qp) g \ > |Im(gp)J: 



(?) g (?) g >\m e \ 2 >\m,P}) e \ 2 = (^ 



2 



Note that Heisenberg's uncertainty principle is usually proved only for a single 
state g = g in the form of the Robertson inequality Ro > fi 2 S ( ^" 1 /4 for the vari- 
ances Ro and So of the dynamically conjugate variables described by the canonical 
operators qi and p k in TL which satisfy the exact canonical commutation relations 

[qiAk] = 0, [qi,Pk] = i^jfeJ, [pi,p k ] = o. 

A more precise matrix multidimensional generalization of the Robertson inequality 
in terms of the covariances of estimates of an arbitrary family of noncommuting 
operators is proposed in [J| . Note that Robertson inequality implies the uncertainty 
relation 

Re > fr 2 S 1 /4, S = [TrsiSkg ] = S g 

for the unbiased measurements of the unknown expectations 9i = {<ji) g in the canon- 
ical states (13.311 with J Xill (dA) = §i, where qi satisfy the canonical commutation 
relations with the canonically conjugated §k = pk- In this case the unbiasness 

Me [9] = (q) g = Tvg q (9) = (q) Q + 9 = 9, 
simply means that (<?i) = for the state g as 

q{9) =e~ ie »° k qe WkSk =q + 9 
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Every such unbiased measurement has the variance Re > Qe, and among such 
measurements there is an optimal one corresponding to Rg — Qe. It is realized 
by the direct measurement of all qi described by the orthogonal spectral measure 
II (dA) = E (dA) of the commutative family q~i = f XiE (dA). Note that in this case 
p(0) = p, and both Sg = So and Qg = Qq do not depend on 9 in any state g — g g . 

Our analysis extends the Heisenberg uncertainty principle to any unbiased mea- 
surement satisfying {q) e = 9. Note that without unbiasness the uncertanicy relation 
doesn't hold for such dynamically conjugate variables as polar coordinate described 
by the bounded selfadjoint operator — ttI < q < irl and the discrete angular mo- 
mentum s. In this case one can find a state g (e.g. the eigen state of angular 
momentum for which the uncertancy relation is obviously not true as So = and 
Qo < 7r 2 ). There is no good operator q in Tt satisfying the unbiasness condition 

(<lig "■ 

We now consider the general case of the non-commuting generators x k in the 
canonical family H3.ll) . Differentiating (|3.1|1 with respect to p k and comparing the 
result with (|2.5(l . we obtain 

(3.6) h k = e-^ x 4fcX _1 e^ = x k (fl) - Mfel 

where x k = e - ^ Xk Xk , and Mfc = Mfc ~P) 1S tne expectation value of 

x k (P) at the state g — g (/3, /3) : 

H h = Tri fc (p) g (p, P) = 4^ ln X (P, P) ■ 

dp 

The right Fisher information matrix (|2.6|l is therefore the matrix of the covariances 

(3.7) h ik = Tr (xt (P) - Mi) (x k CP) -Vk)*Q (P, P) = 

dp dp 

of the operators x k (P) depending analytically on P (but with not necessarily ana- 
lytic expectations x k at g(p,/3)). The inequality (|2.7() in the neighborhood of the 
point P = can therefore be expressed in the form of the uncertainty relation 

(3.8) R (P, p) > D (p, p) S (P, P)~ X D (p, p) f , 

which establishes an inverse proportionality between the matrix S = [<Sifc (/3, /?)] of 
the covariances 

(3.9) S lk = Trg (P, P) (x l - n t ) (x k - fi k )* 

for the operators x k — x k (0) with the expectations [i k = Trx k g (P,P) and the 
covariance matrix R (/3, P) of the estimates A^ for the functions $i (/9, P) of the 
canonical parameters p k . 

Let us consider the case when the operators x k are the generators of a Lie algebra. 
Suppose the operators x k satisfy a Lie algebra commutation relations 

(3.10) xix k x k Xi = Cj^Xj, 

where C 3 ik are the structure constants. In this case, the operators Xi (ft) in l3.6|) are 
linear combinations of the generators Xi = ii (0) |12|: 

(3.11) x l CP) = K -1 (p)i 
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where K (0) = /3 k C k (tP" c > - l) is an n x n matrix which exists in, at least, 

a certain neighborhood O C C" of the origin j3 = 0, and C k = C\ k are the 
generators of the adjoint matrix representation 

C^Cfe C/^C^ C^Qj 

of the commutation relations (|3.1U|) . Expressing the covariance matrix H of the 
operators i|3.11|l in terms of the covariances l|3.9[) of the generators Xi, we obtain in 
place of 13 .8(1 the exact inequality 

(3.12) R > (DK+S-W) {p,p). 

In the case (|3.3|) . the family g e is unitarily homogeneous with respect to the Lie 
group having Hermitian generators Xf. = Sf. = x* k and canonical parameters 9i. As 
in the case of ()3.5|l . we obtain a generalized uncertainty relation 

(3.13) R 9 > ^KjSo%, 

where K = i8 k C k (e ie » c " - 1) 1 and C fe = fi _1 C fe . In the domain 9 C W l of 
convergence of the series 

m— 1 

the inequality (|3.13|) determines the lower bound of the mean quadratic error of 
measurement of the canonical parameters for the unitary representation e l6kS of 
the Lie group generated by the selfadjoint s k — h~ 1 §k- 

4. Efficient Measurements and Quasimeasurements 

1. In classical statistics, estimates whose covariance matrix attains the minimal 
value, transforming the Cramer-Rao inequality locally or globally into an equality, 
are said to be efficient (locally or globally, respectively). In quantum statistics, 
because of the nonunique generalization of the Cramer-Rao inequality, the concept 
of efficiency, introduced by analogy with the classical concept, loses its universality, 
and the definitions of locally efficient estimates pQ, based on the different 
variants of this generalization are not equivalent. Therefore, we shall distinguish 
efficient measurements (or estimates), for which the invariant Helstrom's bound 
(|2.4(l is attained, from efficient measurements corresponding to the right bound (|2.7J) . 
calling the former Helstrom efficient and the latter right efficient. As we shall 
show here, the concept of right efficiency is more universal: Measurements that are 
globally Helstrom efficient are also right efficient, but not vice versa. We show first 
that Helstrom efficient estimates exist globally for the canonical families of density 
operators l|3.1|l if the operators x k are commuting self-adjoint operators Sfe, and the 
estimated parameters $ (7) are taken to be their expectations 

(4.1) ^( 7 )=Trs^(7)=// fe (7) 

as the derivatives = dlnx/dl k f° r the moment generating function x (7) = 
Tig e yhsk of the canonical states 

(4.2) Qh) = X- 1 (l)e-' kSk/2 Q ei ks »/ 2 
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corresponding to zero imaginary parts Im (3 k = in with \ (/?, /?) = \ {fi + ft) . 
Differentiating the operator-function 1(4.2(1 we find the symmetrized logarithmic 
derivatives g k = s k —/i k with respect to j k . Thus, the symmetric Fisher information 
(12.21) in this case is the matrix of covariances 

(4-3) S ik = Trg (7) (s; - Mi ) (s k - Mfc ) = . 

for the operators s k . However these covariances as the second derivatives of In \ 
are the derivatives d\i i jdy t = dfi k /dj l of (|4.1|) . That defines the matrix D = 
[<9iy<97 fe ] in (21} as 

D( 7 ) = [a Mi (7)/a 7 fe ] =s( 7 ). 

The inequality (|2.4I) therefore takes the form Q (7) > S (7), i.e. [Qi k — Sik] > 0, 
where Q (7) = R (7) is the covariance matrix of the operators q = s realizing the 
unbiased estimates by the joint measurement of §i. One can take the spectral QPM 
II of the family Sj = J xjl (dx) and define these estimates as spectral values x k for 
s k . The covariance matrix R (7) of such estimates obviously achieves its minimal 
value 

R = M,j [(A, - ft) (A fe - fc )] = M fl [(xj - Ml ) (xfc - /i fe )] = S. 
Thus, for the canonical families l|4.2|l with commuting self-adjoint s k there exists a 
Helstrom-efficient estimation A = x of the functions (|4.1|l defined by the canonical 
parameters n k , and this is realized by an a simultaneous measurement of the com- 
muting observables s k . The domain of this efficiency obviously coincides with the 
domain O C R" in which \ (7) < 00 is twice differentiable. It can be shown that 
the opposite assertion holds in the following sense. 

Suppose that the estimates Ai (i.e., the results of a measurement) have, in a 
certain domain, differentiable mean values (a) and the covariances Ri k (a), and 
suppose the matrices R = [Rik (a)] and D = \d , di/da k ] satisfy the following regu- 
larity conditions 

(which are trivial in one-dimensional case). Then one can introduce the canonical 
parameters j k by setting j k (oto) = for an ao at which $ («o) = 0. 

It is readily verified that for a family of density operators g (7) of the canonical 
form Q4.2[l the regularity conditions are satisfied for the efficient measurement of 
*& k = /i k (7) as in this case 

R( 7 ) = S( 7 ), D( 7 )=S(7) 

and therefore (R _1 D) = I. The proof of the opposite assertion, that if the regular- 
ity conditions are satisfied, global Helstrom efficiency holds only for the canonical 
families 14.2(1 . is given in the Appendix for the more general complex situation. 
Thus, 

Theorem 1. Under the above regularity condition the inequality \2.$ becomes an 
equality in the domain O C R" iff the density operators g (a) have the canonical 
form \4-2\ ), where s k , k — l,...,rt, are Hermitian commuting operators in TL, 
and the canonical coordinates 7 are functions of the parameters a defined by the 
equations 

d 

—glnxiv =$k(a), k = 1, . . . ,rc. 
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The optimal estimation in this case reduces to the measurement of the Hermitian 
operators §k described by their joint spectral resolution of identity, and the minimal 
mean square error is determined by the matrix of their covariances H-'Jty . 

2. Suppose that in a domain O C C™ of some complex coordinates a — (a fc ) 
the unbiased estimates Xk have mathematical expectations (a) and covariances 
Rik {ex, a) satisfying the regularity conditions 14.4(1 

(4-5) ^(R-W ^-^(R-W , _^R- 1 D = 0. 

y ' da lK > k da k v h ' da k 

(which simply means in one-dimensional case the analyticity dR~ 1 D/da = 0). 
Then, as in the real case, one can introduce the canonically conjugate parame- 
ters f3 k = (3 k (a) as analytic functions satisfying the the equations 

dftldct = (R _1 D)* , (3 k (a ) = 0. 

and the functions f3 k (a) are analytic by virtue of the condition (|4.5I) . 

Theorem 2. Under the above formulated regularity conditions, the inequality \2. 61) 
becomes an equality if and only if the family {g(a, a) ,a £ O} has the canonical 
form i(ff.i() . where g — g (0, 0), the operators Xk, k = 1, . . . , n, have simultaneously 
in TL the right eigen QPM 

(4.6) i=/n(dx), x fc n (dx) = K k n (d^) , x = (x 1 ,...x n )ec n , 

and the canonical parameters j3 k , k = 1, . . . , n are defined by the equations 

4.7 = ^ (a, a), a e O. 

8(3 

The optimal estimation in this case reduces to a quasimeasurement of the non- 
Hermitian operators Xk described by the resolution of the identity \4-.b\j , and the 
minimal mean square error is determined by the matrix of the covariances 

(4.8) H ik =Tr e {xi-fi i ){x h -■&*)*. 

The sufficiency is proved as in Section |21 Using the invariance of the right 
bound 1(2. 7JI under the analytic transformations a i— > /3, we choose as the displaced 
a k determining this bound the canonical parameters (3 of the family of density 
operators 1(3. The elements d-di/dj3 k of the matrix D with allowance for d u = 
dhix/djf then coincide with the elements (I3.7|l of the matrix H. Since the operators 
Xk commute in accordance with (14.6(1 , 



*A = y wiCd*)^, 

we have "dk = M/ci H — S, where [i k are the mathematical expectations of Xk and 
S is the covariance matrix ((3.91) of these operators. Therefore, the inequality ((2.7(1 
takes the form R > S. It remains to show that the measurement described by 
the resolution of the identity ((4.6(1 leads to an estimation for which R = S even in 
the case when the operators Xk do not commute with their Hermitian conjugates: 
x*Xk x* k Xi (which is the case for a nonorthogonal resolution 1(4.6(1 1. For this, it is 
sufficient to take into account the representation 

(4.9) x l — J \tiU(dx), XiX* k = /xiXfell(dx), 



14 



V P BELAVKIN 



obtained by integrating the equations in (|4.6(l x E C" and also the conjugate equa- 
tion II (d>r) x£ = iir/cll (dxr). Because of ()4.9|l . the covariances 



(4.10) i? ifc = / fa - 0<) (x fe - tf fc ) Trgn (dx) 

of the estimates A& = obtained on the basis of the quasimeasurement of the 
operators q~k = Xk coincide with the covariance Hik of these operators, which proves 
that this generalized measurement is efficient for the density operators 1|3.1|) . The 
proof of the opposite assertions of Theorem 2 follows from the very derivation of 
the inequality (|2.7|) and is given in the Appendix. 

3. Thus, the condition of (right) efficiency requires the existence of commuting 
operators that have a joint right spectral resolution and play the role of sufficient 
statistics, which we call right-efficient. At the same time, it is sufficient to restrict 
the study of these operators to the minimal subspace generated by the domains 
g ((3,(3) H with density operators g ((3,(3) for (3 G O. further, if one considers only 
real values of the parameters $k ((3, J3) , the optimal estimation can be described by 
non-Hcrmitian and noncommuting (with the conjugate) operators of right-efficient 
statistics and is not therefore Helstrom efficient. However, estimates that are Hel- 
strom efficient correspond, in accordance with Theorem ^ to the special case of 
right efficiency for which the operators Xk are Hermitian. If the operators Xk in 

1(1 are not Hermitian but commute with the Hermitian conjugates, the right ef- 
ficient estimates also coincide with the complexified Helstrom efficient estimates 
However, the commutativity XkX* = x*Xk need not hold. 

Example. Suppose Xk — p>k (^)j where ip k are entire functions C — > C, a = 
(a,, . . . ,a r ) are boson annihilation operators satisfying the commutation relations 

[a i: aj] = 0, [a J: a*] = 

It is well known that the operators a have right eigenvectors |a) € H, a G C r , that 
define a nonorthogonal resolution of the identity: 



f r 1 

1= / \a) (a\ — dReaidlmai, \a) = a, |n 



Obviously, the operators x, — cp (a) also have a right eigen resolution of the identity 
(|4.5|l . where 

n (d>r) — / S (dx , ip (a)) \a) (a\ JJ^ — dRe aidlmai 

i—l 

(5 (dx, A) is the Dirac delta measure of unite mass at the point A). Therefore, the 

— k 

optimal estimation of the parameters i?fe = d In \l 9(3 of the density operators (|3.1() 
for x — if (a) is right efficient and reduces to a coherent measurement and extension 
of the estimate i? = (p (a) with respect to the result a. For the special case when 
the function (p (a) is linear and the state g is Gaussian, this fact was established 
in 0. 

Note that besides right and left lower bounds one can also consider other, com- 
bined bounds by means of the factorization ^5] d = i?+ +1?-, defining right deriva- 
tives with respect to and left derivatives with respect to An interesting 
question is this: Is the class of efficient statistics exhausted by statistics for which 
at least one such bound can be attained? 
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4. In conclusion, let us consider the question of the (right) efficiency of the 
estimation of the parameters f3 k themselves of the canonical families <|3.1|) . The 
inequality (|2.7I) corresponding to this case d k = P k has the form R > H" 1 where H 
is the matrix of the derivatives (|3.7|1 . Without loss of generality, we shall assume 
that TrxkPQ = 0. 

Theorem 3. The inequality R > H _1 becomes an equality if and only if the 
operators x k in $3.1\) have a right joint spectral measure j4-5\) , the generating 
function of the moments y3. £)) of these operators in the state g is Gaussian: 

X (/?,/?) = exp ^J3' Hikf3 k ^ , where does not depend on (3 and 0, and the un- 
biased estimates X k — X k (x) are taken to be linear functions X k = H kl Xi of the 
results x k of simultaneous quasimeasurement of the observables Xk ■ 

The proof of the sufficiency of these conditions for the existence of the right 
efficient estimation is obvious: From the fact that the matrix H coincides with 
the covariance matrix S of the operators Xk it follows that the covariance matrix 
R = H-iHI-r 1 is equal to H" 1 . 

The necessity follows from the necessary conditions of right efficiency of Theorem 
13 according to which the family g (/?, p) must also have the form 

(4.11a) Q {f3,0)=X- 1 e^ h "Q o e^\ 

where x (/?,/?) = TtqqqP* x ef** x * , -^-\n\ = (3 k ', and the operators x k have the 
joint right resolution of the identity 

i = J n (dx) , x k u (x) = x k n (dx) , x = (x k ) e C". 

Comparing l|4.2[l and l|3.1|l . we obtain ft k x k = j3 k x k , whence 

h = H ki (3\ X (/?, 0) = p l H lk p k , x k = H kl Xl . 
Theorem 3 has been proved. 

5. Appendix 

1. Let us proof the inequality l|2.7|l . First consider the one-dimensional case. 
Let q be an operator in TL for which 

(5.1) Trqg (a, a) = d (a, a) . 

Differentiating l|5.1|) with respect to a and using the definition (|2.5|l and the nor- 
malization condition Trg (a, a) — 1, due to which Trgh* = 0, we obtain 

d-d 

= Trg(q-d)h*. 

da 

Since the covariance Tig (q — h* satisfies the Schwarz inequality 

(5.2) T \rg (q - i?) h*\ 2 < Tr [g (q - &) (q - #)*] Tr \ghh 

which is the condition of non-negativity of the determinant of the 2x2 matrix of 
covariances Tigh^l, i = 0, 1, where ho = (q — $) , h\ = h, we can write 

2 / 



(5.3) Trg {q - 1?) (q* - ti) > 



d<d 
da 



Tr ghh* , 
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This inequality obviously establishes a lower bound for the variance of the estima- 
tion of the parameter = ■& (a, a) in the class of ordinary measurements described 
by normal operators q. However since the normality condition qq* — q*q was not 
used in the derivation of (|5.3fl . this bound gives a lower bound for the variance 
of any unbiased estimation of Indeed, if II (dA) , A E C is a QPM describing 
the unbiased estimation as a generalized measurement in Ti., then the operator 
q = J All (dA) satisfy the condition 1(5.1(1 . From the Hermitian positivity 

(5.4) (A-g)II(dA) (A-g)* > (II > 0) 
it follows that / |A| 2 II (dA) > qq* , and 

(5.5) y|A-g| 2 n(dA)>(A- ( 7)(A- ( 7)*. 

Taking the mathematical expectation of both sides of (|5.4I) and bearing in mind 



that the variance R = M$ 



of the estimation $ is 



iA-*r 

R = Trg J |A - d\ 2 n (dA) , 



we obtain in conjunction with ((5.3(1 

(5.6) R>Tr g (q-#)(q-#y > \d\ 2 /g, 

where we have denoted d = d$/da, g = Tvghh*. Thus, for the one-dimensional 
case the inequality 1(2.7(1 has been proved. 

2. Equality can be attained in 1)5.5(1 if, first, the expectations of the two sides of 
((5.5(1 coincide and, second, the Schwarz inequality becomes an equality. The first 
condition actually establishes equality in 1(5.4(1 . More precisely: 

Lemma 1. Suppose the ranges g(a,a)TL of density operators {g(a,a) : a E 0} 
generate the whole ofH. Then the equality TvgR = for any non-negative definite 
operator R in Ti and all a E O implies that R = 0. 

It is sufficient to show that in H there is no vector \x) of the form |x) = g 1 / 2 
for which (\\ R \x) 0- But this follows from the inequality 

Trg 1 / 2 Rg 1 / 2 > (</>| g 1 ' 2 Rg 1 ' 2 |^> . 

which holds for any non-negative R when (ip \ ip) = 1. 

Applying this result to the operator R equal to the difference of the right- and 
left-hand sides of ((5.5(1 . we find, under the conditions of the lemma, that equality 
holds in ((5.5|) only if 

(A - q) n (dA) (A - qf = 0, or qll (dA) = An (dA) . 

This proves that for the existence of right efficient unbiased estimation in some 
domain O 3 a it is necessary to have an operator q with a right-eigen QPM in the 
subspace generated by the subspaces g (a, a) H, with Trqg (a, a) = A. In the case 
of real spectrum A E K such an operator q is obviously selfadjoint. 

The second condition of equality in ((5.6|) is equivalent to the condition of linear 
dependence g (q — A) = tgh, where t = d/g is a constant. Setting 

ts = q - & (0) 

we obtain the equations 

dg/da = g{s — //) , dg/da = (s — ji) g 
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where t[i = $ (a)— ?9 (0). Its solution of these equations with the boundary condition 
£>(0,0) = g has the canonical form The the operator q = ts + i? (0) should 

have right-eigen QPM, so the operator s should. This proves for the one-dimensional 
case, the necessity of the canonicity of the density operators g (a, a) for the existence 
of the right efficient estimation formulated in Theorem 2. In the Hermitian case 
x* = x, this also proves the necessity of Theorem 1. 

3. A multidimensional generalization is obtained from the one-dimensional case 
by taking 

q - A = (qi - Xi) rf, h = h k (, , 
where rf , i — 1, . . . ,m, a k , k = 1, . . . , n, are arbitrary complex numbers. Remem- 
bering that then 

Tig (q - A) h* = rf ^-r£ fc , 
dor 

we obtain from (|5.2|l for £ = (H _1 Dt) . rf the second of the inequalities 

RikvV > Trg (qt - Xi) (q k - A fc )* v'v" > (DH _1 D t ) ifc ffrf, 
which holds for arbitrary qi for which Tr gq*i = i%. Setting 

qi = J X t U (dA) , where / n (dA) = i, A £ C m , 

is the resolution of the identity describing the estimator Xi — Ki , and applying the 
inequality l|5.5|l for q = qifj 1 , X = Xiff , we obtain for the matrix R of the covarianccs 
of $i satisfying the first of the inequalities H5.6JI . whence l|2.7|l follows because rf is 
arbitrary. 

The inequality Ij2.6|l becomes an equality for a €E O only if 

q t n (dA) = A. t n (dA) , and g (q, - A,) = f k dg/da k , where t\ = (DH -1 )^, , 

whence with allowance for T = [T^] to be constants nondegenerated matrix we 
obtain (pTT|) . 
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